Lepton Number-Driven Sterile Neutrino Production in the Early Universe 
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We examine medium- enhanced, neutrino scattering- induced decoherent production of dark matter 
candidate sterile neutrinos in the early universe. In cases with a significant net lepton number we 
find two resonances, where the effective in-medium mixing angles are large. We calculate the lepton 
number depletion-driven evolution of these resonances. We describe the dependence of this evolution 
on lepton numbers, sterile neutrino rest mass, and the active-sterile vacuum mixing angle. We find 
that this resonance evolution can result in relic sterile neutrino energy spectra with a generic form 
which is sharply peaked in energy. We compare our complete quantum kinetic equation treatment 
with the widely-used quantum Zeno ansatz. 
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I. INTRODUCTION 

In this paper we examine the physics which determines 
the relic densities and fossil energy spectra of neutrinos 
in the very early universe. The existence of new kinds of 
electroweak singlet "sterile" neutrinos, Vg, with subweak 
interactions remains a controversial possibility. The ex- 
perimental establishment of nonzero neutrino rest masses 
fuels legitimate speculation on this issue. 

Moreover, this issue could be important because ster- 
ile neutrinos with rest masses ~ 1 keV are interesting 
dark matter candidates. The minimalist model for this 
dark matter scenario is where sterile neutrinos are cou- 
pled to active species via mass terms that give vacuum 
mixing. In turn, the simplest model along these lines 
is where active neutrino scattering-induced decoherence 
gives rise to a relic density of sterile neutrinos [l| . 

These simple models, however, are challenged. Vac- 
uum active-sterile neutrino mixing, which is necessary 
to produce the sterile neutrinos from an initially purely 
active neutrino population, also enables a non-GIM sup- 
pressed radiative decay channel [vg — > Ve,p,,T+"t) for ster- 
ile neutrinos. As a consequence, the best constraints and 
probes of this sector of particle physics come from the 
modern x-ray observatories 0, d, |j, H, 0, 0, H, H, [13, EH ■ 
These constraints arguably eliminate the simplest models 
for decoherence-produced sterile neutrino dark matter. 

However, there remain processes for sterile neutrino 
production which can produce the correct relic density for 
dark matter yet evade all current bounds. These viable 
models include, for example, sterile neutrino production 
associated with inflation [ij] , Higgs decay [ij, jlal , 
and lepton number-driven medium enhancement [3|, |l6| . 

The concept of lepton number is a slippery one if 
sterile neutrinos exist and mix in vacuum with ac- 
tive species. In essence, lepton number can be cre- 
ated [III or destroyed [l3, [13l by active-sterile neu- 
trino interconversion. The current observational bounds 
^9, 21, 22, 23, 24, 25, 2a, H] on the electron, muon, and 
tau lepton numbers of the universe are relatively poor, 
at least on the scale of these quantities required to affect 



active-sterile neutrino conversion. Therefore, it may be 
important to examine in detail how lepton number influ- 
ences scattering-induced decoherent production of sterile 
neutrinos in the early universe. This is the objective of 
this paper. 

In what follows we investigate how the coupled inter- 
play of sterile neutrino production, lepton number de- 
pletion, and the expansion of the universe determine the 
relic sterile neutrino densities and energy spectra. In 
Sec. In] we describe in-medium active-sterile resonances 
and how these influence decoherent sterile neutrino pro- 
duction. In this section we also discuss the relationship 
between the quantum kinetic equations and the quantum 
Zeno approximation. In Sec. IIIII we discuss the approach 
we employ to solve for lepton number and Vg spectral 
evolution. In Sec. IIVI we discuss results. Conclusions are 
given in Sec. IVl 



II. STERILE NEUTRINO PRODUCTION IN 
THE EARLY UNIVERSE 



There is a long history to the general problem of 
the production of sterile neutrinos i^s from an initial 
population of purel y a ctive neutrinos Va (a = e,fj,,T) 
[ll,H, [H, [13, [H, [29ll30j . In low density environments in 
supernovae and the early universe, the dominant Vg pro- 
duction channel will be the coherent Mikheyev-Smirnov- 
Wolfenstein (MSW) process [3l|, H^l- However, neutrino 
propagation may not be coherent in the higher density 
regions of core collapse supernovae and the pre-neutrino- 
decoupling epoch in the early universe. In these envi- 
ronments, scattering-induced decoherence is the princi- 
pal way in which sterile neutrinos are produced. Here we 
concentrate on the decoherent production channel since 
this will dominate the production of dark matter candi- 
date sterile neutrinos. 
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A. The Quantum Kinetic Equations 

The evolution of a system of active and sterile neutri- 
nos can be formulated in terms of the density operator 
for a neutrino with given scaled momentum e = p^/T^ 
with the neutrino momentum and T the temperature 
in the early universe. (Here we use natural units in which 
Ti = c — — \.) The density operator for scaled mo- 
mentum e is 



(1) 



where the summation is performed over neutrino flavor 
states, e.g., \vi) with i = e, ^,t, s and py (e, t) are the cor- 
responding density operator matrix elements which are 
momentum and time dependent. Here we take the tem- 
perature T(t) to be a function of the Friedman-Lemaitre- 
Robertson- Walker time coordinate t, defined by the so- 
lution of the zero-curvature Friedman equation 



87rp(t) 
3mpL 



(2) 



where a{t) is the scale factor, p{t) is the total energy 
density of the universe, and mpL is the Planck mass. 

For illustrative purposes, we consider 2x2 vacuum 
mixing between an active and sterile neutrino with a one 
parameter (vacuum mixing angle 6) unitary transforma- 
tion between the weak eigenstates |i^q) (a = e,/i,T), 
and the energy-mass eigenstates 1^2) ■ 



\Va) = cosflt^i; + smfc'|t/2/; 
\vs) = - sin^li^i) -I- cos6'|j/2) 



(3) 



In this 2x2 formalism, it is convenient to decompose the 
density operator as 

p{e,t)^]^P,{e,t)[l + V{e,t)-cT], (4) 

where cr is the Pauli spin operator, Po(e, t) is a normaliza- 
tion factor proportional to the total number of neutrinos 
(active and sterile) with scaled momentum e, and P(e,t) 
acts as a polarization vector in weak isospin space. 
The diagonal matrix elements of the density operator, 

/„(e,t) = p„„(e,i) = ]^P^{t,t)[\ + P,{e,t)] (5) 

/,(e,t) = p,,(e,i) = ipo(e,t)[l-P,(e,t)], (6) 

are proportional to the number density distributions of 
the neutrino species. The density operators are normal- 
ized so that the neutrino number densities are given as 
products of the above matrix elements and the zero chem- 
ical potential Fermi-Dirac distribution: 



27r2 e<^ + 1 ' 
27r2 + 1 ' 



(7) 
(8) 



The quantum kinetic equations can be derived [331 134| | 
from the time evolution of the density operator. In a 
homogeneous and isotropic universe, these equations can 
be given in terms of the time evolution of the functions 
Po(e,^) and P(e,i) as defined in Eq. (g]), 



V(e,i)xP(e,i) 



^lnP„(.,0 



d 



Die,t) + —\nPoie,t) 



d_ 

m 



Poie.t) = R{e,t), 



P^(6,t) (9) 

(10) 



where P ± — PxX + Pyij , and the functions V, £>, and R 
are defined below. 

The vector V(e, t) corresponds to the coherent quan- 
tum mechanical evolution of the neutrino states in the 
early universe, and is given by 



Y{e,t) 



2eT{t) 



where Sm? — 



(sin26'£ - cos26'z) + Va(e,t)z, (11) 

— ml is the difference of the squares 
of the vacuum neutrino mass eigenvalues. The first term 
represents vacuum neutrino oscillations, while the sec- 
ond term introduces matter effects through the forward 
scattering potential [sH, [sHi 

K(e,i) = ^^i^^^ GpCa^m^it) - r^GleT'it), (12) 

where ({3) sa 1.20206, Gp is the Fermi constant, and 
is a dimensionless coefficient which depends on the 
number of charged lepton degrees of freedom and can be 
found in Refs. 0, [l^. The potential lepton number is 



Ca{t) = 2L^{t)+yLp{t) 



The individual lepton numbers are 



Lait) 



At) 



(13) 



(14) 



where n^^, np^, and are the neutrino, antineutrino 
and photon number densities, respectively. Here we ne- 
glect contributions to the forward scattering potential 
from neutrino-baryon and neutrino-charged-lepton inter- 
actions since we will consider lepton numbers much larger 
than the baryon to photon ratio H, [l^ . 

The decoherence function D(e, t) and the repopulation 
function R{e,t) can be simplified by assuming thermal 



equilibrium of the background plasma [37[. The deco- 
herence function corresponds to the loss of coherence re- 
sulting from collisions with particles in the early universe, 
and is proportional to the total scattering rate Ta{e, t) of 
neutrinos Vn, 



Die,t) = -r„(e,t). 



(15) 
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Using the assumption of thermal equilibrium, the total 
scattering rate can be written as 



(16) 



where we have neglected corrections of order the lep- 
ton number La. The numerical coefficient ya{t) pri- 
marily depends on the number of relativistic particles 
with weak charge that are populated in the thermal seas 
of the early universe at epoch t. For example, at tem- 
peratures 1 MeV < T < 20 MeV the total scattering rate 
mostly stems from interactions with other neutrinos and 
pairs, so that j/e ~ 1-27 and j/^^,- « 0.92 However, 
at higher temperatures there can be appreciable popu- 
lations of other charged leptons and quarks which will 
increase the total scattering rate. 

The repopulation function R{e, t) dictates the evolu- 
tion of Po{e,t). Since Pq is proportional to the total 
number of neutrinos with scaled momentum e, the re- 
population function corresponds to scattering into and 
out of neutrino states {i>a and Vs) with scaled momen- 
tum e. Assuming that each populated species, with the 
exception of i^a and Vs, has a thermal spectrum, the re- 
population function can be written as 



i?(e,t) =r„(e,t) 



+ 1 



(17) 



where rja is the degeneracy parameter (ratio of the chem- 
ical potential to the temperature) associated with the 
Fermi-Dirac spectra of i>a and Da with lepton number 



La — 



12C(3) 



(18) 



This form of the repopulation function is valid for tem- 
peratures above the neutrino decoupling temperature, 
where the active neutrinos are able to efficiently exchange 
energy and momentum with the plasma of the early uni- 
verse. 

The effect of the repopulation function is to drive the 
distribution of the active neutrino toward a Fermi-Dirac 
spectrum consistent with a lepton number La, 



na{e,t) cx fa{e,t)- 



(19) 



B. The Quantum Zeno Ansatz 

The quantum kinetic equations, even in their simplified 
form with the assumptions of homogeneity and isotropy 
of the early universe and thermal equilibrium in the back- 
ground plasma, comprise a system of coupled nonlinear 
integro-partial differential equations that are difhcult to 
solve. Past works have employed the quantum Zeno ap- 
proximation to circumvent these difficulties [3, [1^ . We 
write the quantum Zeno ansatz as a Boltzmann-like ki- 
netic equation. 



where the effective matter mixing angle 9m is defined by 



sin" 20M(e,i) 



|V(6,i)P 



sin^ 261 



and the neutrino oscillation length is 

£„(e,t) = |V(e,i)ri 



\2eT{t) 



sin 29 



ra{e,t)em{e,t) 



[fa{e,t)-f4e,t)], 



( 5m? 



(20) 



9 / 2eT(t) , , 
sin^26l-f ( cos 26* — V^y^ e,i 



(21) 



\2tT{t) 



cos29 -Vaie,t) 



-1/2 



(22) 



The de/dt term takes into account the change in the 
scaled momentum of a neutrino propagating along its 
world line through the expanding early universe. Using 
the conservation of comoving entropy and assuming ra- 
diation domination, this is 



de 
dt 



_^dg 
3(7 dt ' 



(23) 



where g is the effective statistical weight in relativistic 
particles and is calculated with a weighted sum of the 
bosonic (gb) and fermionic (gy) degrees of freedom. 



E 



9b 



■ E 

fermioiis 



9f- 



(24) 
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For example, at temperatures 1 MeV ^ T < 20 MeV the 
plasma of the early universe consists of photons {gi, — 2), 
pairs {gf — 4), and neutrinos and antineutrinos of 
all three active flavors (<?/ = 6), so g ^ 10.75, where we 
neglect corrections of order the lepton numbers. 

Using the quantum Zeno ansatz instead of the full 
quantum kinetic equations eases computational demands 
by allowing calculations to be performed with distribu- 
tion functions (amplitude squared) instead of with quan- 
tum amplitudes. The former approach avoids the com- 
putational pitfalls of rapidly varying complex phases that 
plague the latter. In Sec. IIVI we will compare the ster- 
ile neutrino production rates calculated by the quantum 
kinetic equations to those implied by the quantum Zeno 
ansatz. 

Inspection of the quantum Zeno ansatz, Eq. (|20p. sug- 
gests that maximal sterile neutrino production should 
occur when sin^ 26m = 1- This corresponds to the reso- 
nance condition in the coherent MSW process, 

Sm^ cos 26 = ^^P^ Gf/:e,esT^ - 2r„G|eLr'- (25) 

The solutions to Eq. ([^S]) are the resonant scaled mo- 
menta 



.(26) 

Therefore, resonances occur when the following condition 
is met: 

1^(^)1 m > ^^r^5m^ COS26. (27) 

In Sec. IIVI we will discuss the results of our numeri- 
cal calculations and the effects of these MSW-like res- 
onances and this resonance condition. The production 
of sterile neutrinos from an initial active neutrino distri- 
bution reduces the potential lepton number C. Further- 
more, the expansion of the universe results in a decrease 
of the temperature T. These two trends together ensure 
that the resonance condition, Eq. (P7|) . will be violated 
at some point. However, the enhancement in sterile neu- 
trino production associated with this MSW resonance is 
suppressed by the quantum Zeno effect. 

The quantum Zeno effect is the result of scattering- 
induced decoherence interrupting the accumulation of 
quantum phase, suppressing quantum transitions be- 
tween two discrete states [38|. In the quantum Zeno 
ansatz, this effect is represented by the multiplicative fac- 
tor, 



The ratio of the neutrino oscillation length, 1^, to the 
mean scattering length, F"^, determines the level of sup- 
pression stemming from the quantum Zeno effect. 



If the neutrino oscillation length is much larger than 
the scattering length, £„i ^ F"^, quantum coherence is 
lost before significant quantum phase can be accumu- 
lated, suppressing the transition. This is the case at 
the MSW resonances because at resonance the neutrino 
oscillation length is maximal. The suppression of ster- 
ile neutrino production at the MSW resonances means 
that maximal sterile neutrino production occurs at a time 
that does not correspond to an MSW resonance. How- 
ever, since we will be considering small values of sin^ 29, 
the MSW resonances have small widths and as a result 
we nevertheless have maximal sterile neutrino production 
near to these MSW resonances. 

Far from resonance, there is little suppression from the 
quantum Zeno effect. This results from a coincidence in 
the functional form of the forward scattering potential 
and total scattering rate. At temperatures much larger 
than the resonance temperature, the product Taim ap- 
proaches a constant much less than 1, while at tempera- 
tures much lower than the resonance temperature, Ta£m 
rapidly falls as the temperature of the expanding universe 
decreases. 

There is another possibility for maximizing the sterile 
neutrino production rate. This occurs when the poten- 
tial lepton number has decreased and/or the temperature 
has fallen to a point where the resonance condition, Eq. 
(P7)) . is no longer satisfied and there are no MSW reso- 
nances. In this scenario, sin^20M reaches a local maxi- 
mum but is less than 1, which means that the neutrino 
oscillation length tends to be much smaller than in the 
resonant case. Additionally, this occurs at later times, 
and thus lower temperatures, during the expansion of 
the universe, which means the product Ta£m is signif- 
icantly smaller than in the resonant case. As a result, 
the overall sterile neutrino production rate can be higher 
than in the resonant case and would have a broader peak 
than in the resonant case. With the right conditions, this 
nonresonant production of sterile neutrinos can become 
significant. 

Throughout this section we have discussed the trans- 
formation properties of ^ i^s using the quantum ki- 
netic equations and the quantum Zeno ansatz. However, 
to consider the general problem of neutrino transforma- 
tion in the early universe, we must also include its CP- 
counterpart, Da ^ Dg. The quantum kinetic equations 
and quantum Zeno ansatz are similar for the j^q, Vs sys- 
tem with a few alterations. The forward scattering poten- 
tial Va [Eq. p2)) ] is adjusted by replacing Ca with — £„, 
and the total scattering rate F^ [Eq. pB]) ] is replaced by 
one appropriate for antineutrinos, Fq,. However in the 
early universe, to lowest order, these two scattering rates 
are equal, Fq, ~ r„. 

In this paper we will consider large lepton numbers 
which helps to avoid the necessity of considering both 
neutrinos and antineutrinos. For a large positive lepton 
number, all of the action occurs in the neutrino sector. 
The effective matter mixing angle for antineutrinos is sig- 
nificantly less than their vacuum mixing angles, while 
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neutrinos may experience the MSW-hke resonances dis- 
cussed above. For large negative lepton numbers, the 
same is true, but for the CP-counterpart. 



III. CALCULATIONS 

Calculating the evolution of the active and sterile neu- 
trino distribution functions remains a daunting task. In 
addition to the quantum Zeno ansatz, Eq. (j20[) . we must 
elucidate the lepton number evolution. We assume that 
the active neutrinos are in thermal equilibrium in the 
early universe and that active-sterile neutrino conversion 
is the only flavor changing interaction. The former as- 
sumption of thermal equilibrium is a good one before 
weak decoupling where the high temperature of the uni- 
verse (T > 3 MeV) ensures that the weak interaction time 
scale is much shorter than a Hubble time, the time scale 
for the expansion of the universe. The latter assump- 
tion amounts to assuming that there is no exotic physics 
other than active-sterile neutrino oscillations. For the 
epoch before weak decoupling, the time evolution of the 
lepton number Lq, is 



4C(3)yo "^'e^ 



d 



1 dt 



(28) 



Together, the quantum Zeno ansatz [Eq. (PO)) ] and lep- 
ton number evolution [Eq. (j28[) ] form a formidable system 
of integro-partial differential equations. To simplify the 
situation, we assume that the effective statistical weight, 
g, is constant. Strictly speaking, this is not the case in the 
early universe because species fall out of thermal equilib- 
rium as time progresses and the temperature decreases. 
However, a well accepted detailed history g{t) does not 
currently exist. As a result, for illustrative purposes, we 
will assume that g is constant. 

By making this assumption we benefit in a number 
of ways. The temperature evolution of the universe, 
T{t), is greatly simphfied. As different species fall out 
of thermal equilibrium, particle-antiparticle pair annihi- 
lation reheats the universe, resulting in a complicated 
temperature evolution. The temperature evolution, as 
given by the Friedman equation, the conservation of co- 
moving entropy, and assuming radiation domination, is 



1 dT 
T~dt 



47r3 

Is" 



1/2 



9 



1/2 



T2 



mpL 



(29) 



Since T(t) is a monotonically decreasing function of time, 
we can invert it to find t{T). We are then able to perform 
our calculations as a function of temperature. 

In addition, a constant g means that, by Eq. (|23p . 
de/dt — Q. As a result, the quantum Zeno ansatz be- 
comes an ordinary differential equation for a family of 
sterile neutrino distribution functions, parameterized by 
e. The only coupling between the different values of e 
that remains is found in the lepton number evolution, 
Eq. 



We choose a value of g that is representative 
of the early universe before the QCD transition, 
170MeV < T < IGeV. In this epoch, the plasma of 
the early universe consists of photons (gb — 2), giuons 
(f^f, = 16), and /i^ pairs [gj — 8), neutrinos and an- 
tineutrinos of all three active flavors {gf = 6), and up, 
down, and strange quarks and antiquarks {gf = 36), so 
g w 61.75, where we neglect corrections of order the lep- 
ton numbers and baryon to photon ratio. At higher tem- 
peratures the other quarks come into the picture, but 
at this point the time scale for expansion in the early 
universe is so fast that there is little net effect from the 
difference in g values. 

We use the forward scattering potential given in Eq. 
([T^ . which is accurate at times after the QCD transition 
and for lepton numbers much larger than the baryon to 
photon ratio. However, we are interested in times be- 
fore the QCD transition where giuons and free quarks 
thermally populate the early universe. The proper form 
of the forward scattering potential remains an open area 
of research [s^ [Z^l . The results presented in this paper 
will remain relevant unless the total forward scattering 
potential arising from the quarks is many orders of mag- 
nitude larger than the forward scattering on color singlet 
baryons and mesons with the same net baryon number. 

On the other hand, the total scattering rate is greatly 
influenced by the additional species thermalized before 
the QCD transition. However, since neutrino-quark in- 
teractions are not well studied, the total scattering rate 
in the QCD epoch remains a question mark. We use the 
values given in Sec. Ill Al although our calculations deal 
with a different epoch in the early universe. In addition, 
we calculate the effects of an increased scattering rate, by 
assuming that this rate is proportional to the number of 
relativistic particles with weak charge populated in the 
early universe. 

We solve the system of integro-differential equations 
posed by the quantum Zeno ansatz for different values of 
e, and the lepton evolution. Using Eq. ([M]) . we change 
variables to follow the sterile neutrino distribution func- 
tions and lepton number as a function of the declining 
temperature of the expanding universe. The initial con- 
ditions are chosen so that there are initially no sterile 
neutrinos, fs = 0, and the initial lepton number is a free 
parameter, but we assume that there is an equivalent lep- 
ton number in each of the active neutrino species. The 
initial (final) temperature for the calculations is chosen 
to be high (low) enough that the choice of this has no ef- 
fect on the final outcome (sterile neutrino spectrum and 
final lepton number). 

For illustrative purposes, we choose to work with elec- 
tron neutrinos as the sole active species which mixes with 
sterile neutrinos {a — e). We do not include mixing 
with or Vr] this includes active-active mixing with 
Ve and among themselves, and active-sterile mixing with 
Us- This is unrealistic, since the neutrino mixing ma- 
trix could be quite complicated. However, our simplistic 
model will at least serve to illustrate trends in sterile neu- 
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FIG. 1: Results for our calculation using the quantum Zeno ansatz. In this calculation, we take the sterile neutrino rest mass 
rus = 64keV, vacuum mixing angle sin^ 26 = 10^^", and initial lepton number Leo = 1-1 x 10~^. The upper left panel shows 
the cumulative sterile neutrino production history in terms of its closure fraction in the current epoch, Qsh^- Below that, in 
the lower left corner, is the lepton number history. The upper right panel displays the final i/s spectrum as a function of scaled 
momentum e. In the lower right panel is a plot that tracks the resonant values of e as our calculation evolves in temperature. 



trino production and lepton number depletion [3, [l9l | . 

In addition to the numerical procedure outlined above, 
we also solve the quantum kinetic equations, Eqs. ([9l- fT0| . 
for a given value of e. Using the lepton evolution calcu- 
lated with the quantum Zeno ansatz, we compare the 
evolution of the sterile neutrino distribution function be- 
tween the two methods. The principal difficulty in solv- 
ing the quantum kinetic equations is the need to resolve 
rapidly varying complex phases. To ease these computa- 
tional demands, we employ an eighth-order Runge-Kutta 
method 41] which allows the use of a larger time step to 
solve the problem. 

IV. RESULTS 

The particulars of phenomena revealed by the calcu- 
lations described in the last section depend on the as- 
sumed values of sterile neutrino mass TOs, vacuum mix- 
ing angle 9, and initial lepton number, L^o, in each 
active neutrino species. (As discussed above, we take 
Leo = Lf^Q — LrQ.) Some ranges of the sterile neu- 
trino parameters can be ruled out by the x-ray obser- 
vations 0, i, H i, i, 0, i, i, [H or by observations 
of smaller-scale large scale structure {e.g., the Lyman-a 
forest) coupled with structure formation calculations 



[i, m, m, Bounds on the lepton numbers from 

big bang nucleosynthesis considerations [U, [H, |23| 
are much larger than those employed in our calculations. 
It is important to note that these lepton number lim- 
its constrain the final lepton numbers, not the initial 
lepton numbers. However, other parameter values re- 
main unconstrained by these considerations and, there- 
fore, are still "in play" for providing a relic sterile neu- 
trino density which could comprise some or all of the 
dark matter. The observed dark matter closure frac- 
tion, ilc, is measured to be il.ch^ = 0.1105 ± 0.0039 where 
h = iifp/lOO kms~^ Mpc~^, with Hq the Hubble param- 
eter [isj. 

A set of parameters that produces the correct relic 
mass density to be the dark matter, but is nevertheless 
ruled out by x-ray observations provides particularly in- 
structive cases. These parameters give lepton number 
depletion and sterile neutrino production histories which 
show a variety of possible behaviors associated with res- 
onant and nonresonant evolution. Scenarios associated 
with these parameters give distinctive features in the fi- 
nal relic sterile neutrino energy spectrum which are di- 
rectly attributable to the way the MSW resonance sweep 
couples with lepton number depletion. 

Results for an example set of these parameters are 
shown in Fig.[TJ For this case we have taken = 64 keV, 
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sin^ 29 = 10-1°, and L^o = 1.1 x lO'^. Using the quan- 
tum Zeno ansatz, our calculations yield a sterile neutrino 
relic density which is consistent with the observed dark 
matter closure fraction. 

The panel in the lower right of Fig. [T] shows the sweep 
history of the two MSW resonances {i.e., as in Eq. [26]) . 
The resonant scaled momentum (cros) of each resonance 
has a characteristic behavior. As the universe expands, 
the temperature T drops, and lepton number is depleted, 
these resonances first diverge and then converge in scaled 
momentum. At a temperature ~ 950 MeV, the tempera- 
ture has fallen and the lepton number has been depleted 
to a point where the resonance condition, Eq. (|27p . can- 
not be satisfied and so subsequently there can be no 
MSW resonances. 

The lepton number history and the cumulative ster- 
ile neutrino production history (as measured by its clo- 
sure fraction in the current epoch flsh^) are shown in 
the lower and upper left-hand panels of Fig. [U respec- 
tively. The temperature where MSW resonances cease to 
exist corresponds to an abrupt shift in the sterile neu- 
trino production mechanism, from resonant to nonreso- 
nant production. Roughly 70 percent of the total relic 
sterile neutrino density is produced resonantly. At tem- 
peratures slightly above those characteristic of the res- 
onance cessation event, there is an epoch of precipitous 
depletion of lepton number and concomitant production 
of sterile neutrinos. 

The sharp decline in lepton number and the subsequent 
loss of resonance are a consequence of the coupling be- 
tween the sweep of the MSW resonances and lepton num- 
ber evolution. As the MSW resonances sweep through 
the Ve distribution, and u^s are converted to i^s's, lepton 
number is lost. In turn, the decreasing lepton number 
accelerates the resonance sweep rate. This positive feed- 
back loop results in the precipitous depletion of lepton 
number. This process ends once MSW resonances cease 
to exist, as seen in our calculations. 

This behavior, coupled with the dominance of resonant 
Us production, results in a characteristic set of peaks in 
the relic sterile neutrino energy spectrum 



1' 



(30) 



where Tf is the final temperature in the calculations. 
These features are seen in the upper right panel of Fig. 
[TJ In this figure, there are two distinct peaks seen at 
e w 0.2 and e w 1.15, and a smaller peak at e w 0.5. The 
three peaks are the result of different physical processes. 
Understanding these processes provides physical insight 
into the mechanism of sterile neutrino production. 

Since the magnitude of the lepton number is small 
throughout the calculation, |Le| <SC 1, we can make the 
approximation that /e(e,T) « 1, neglecting corrections 
of order the lepton number. As a result, Eq. ((20l) can 
be solved analytically if we use the lepton number evo- 
lution LeiT) from our calculations. This will provide an 
analytical backbone upon which to base our analysis. 



Using this approximation, the solution to the quan- 
tum Zeno ansatz yields a relic sterile neutrino energy 
spectrum 



^-(^)^;r:i:T U-^^pj-y^ 7(e,T)drn, (31) 



where 



_ 1 re(e,r)sin^2^M(e,T) 
" AT{T) \ + \Tl{e,Tyi^{t,T) 



(32) 



is the Vf, <-!■ Vs conversion rate per unit temperature in- 
terval. In theory, the initial temperature Tq — > oo, but 
in practice is chosen as described in Sec. IIIII 

For a fixed value of e, 7(5, T) is a sharply peaked func- 
tion of temperature T whenever MSW resonances exist. 
This is the result of this function's proportionality to 
sin^20M, which has sharp peaks when sin^ 26'm = 1 at 
the MSW resonances discussed above. The locations in 
temperature of these resonances are implicitly given by 
the solutions of Eq. (|25p for a fixed value of e. 

However, the widths of the peaks in 7(6, T) do not 
correspond to the widths of the associated MSW res- 
onances. Both siv? 29 m{^,T) and £^„(e,T) are sharply 
peaked functions. The resonance width of both func- 
tions is 5T « H tan 20 (see, e.g. Ref. [l^), where H is 
the density scale height at resonance and is defined as 



n 



1 dVe 

5_ _ 8V2C(3) 
T ^ Ti'^S-m? cos 29 



d log C 
dlogT 



(33) 
-1 



In the neighborhood of each resonance, T^i'^/A ^ 1, so 
near resonance 



1 



sin^ 29m {e,T) 



r(T)r,(e,r) ei,ie,T) 



(34) 



The sharp peaks in both sin^ 26'M(e, T) and £l^{e, T), in 
a sense, cancel each other out. This is because the ratio 
of these functions, 



sin^ 29 M{e,T) fSm^sm29 



2eT 



(35) 



is no longer strongly peaked. 

A measure of the width of the resonances in 7(5, T) 
can be obtained by setting Tli'^/A = 1. Far from res- 
onance, Tlel^/4: < 1, so that 7(e,r) cx eT^sm'^29M, 
which means 7(5, T) has the same resonant shape as 
sin^20A./. When Tlif^jA > 1 closer to the resonance, 
7(e,r) (X e~^T^^^, which means the peaked behavior 
is truncated and the location of the maximum of the 
neutrino conversion rate is displaced from the location 
of the MSW resonance to a slightly lower temperature. 
The change in the forward scattering potential required 
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to satisfy this condition is 5Ve « re/2. The width in 
temperature space corresponding to this potential width 
is 



(5T(e, Trcs) 



SVe 




1 dVe 


-1 




res 


Ve dT 


res 



om"' cos 20 

If we restrict ourselves to values of e for which MSW 
resonances are possible (e < 1.15 for the calculations 
done with the parameters in Fig.[T]), then the sharp peaks 
in 7(e,T) allow us to make the approximation: 



-f{e,T)dT ■ 



,,)'5r(e,T,,s,»), (37) 



where gi is a numerical coefficient of order unity, Trcs.i 
is the temperature of the MSW resonance for scaled mo- 
mentum e, and the sum is performed over every reso- 
nance. This approximation breaks down when two reso- 
nances begin to overlap, a situation we will discuss below. 
The product 



7(e,Ti.os)5T ! 



1 5m^ sin^2e 

2 2erros cos 26* 



(38) 



is equal to one-half of the the dimensionless adiabaticity 
parameter defined in Eq. (38) of Ref. [l^, but with the 
flavor off-diagonal potential set to zero, Ber = 0. Note 
also that the density scale height in this work, Ti, is de- 
fined as the width in temperature of the MSW resonance, 
but the density scale height in Ref. [19j], Ti. — T'^Ti., is 
defined as the physical width of the MSW resonance. 

If the gi in Eq. ip7|) are equal to tt, then this result 
is equivalent to the Landau-Zener approximation. How- 
ever, these numerical coefficients depend on the shape of 
the resonances, in particular on the shape of the forward 
scattering potential Ve- On the other hand, the coherent 
Ve Vs probability is only equal to the Landau-Zener 
approximation when the forward scattering potential is 
linear. It is interesting to see the relationship between the 
coherent Ve ~* J^s conversion through the coherent MSW 
process and scattering-induced decoherent Vs production. 
Reference [l^ used this result without this analysis. 

The low-e peak is the product of two opposing factors: 
a Ve spectrum that increases and a I'e — > i^s conversion 
probability that decreases as e increases. At small values 
of e, the number density of Ve with scaled momentum e 
that are available to be converted into Vg is proportional 
to e^. 

The lower right panel in Fig. [1] shows that at these e 
values there are two resonances, one at a high tempera- 
ture (> 2000 MeV) and another at a lower temperature 
near the resonance cessation event (^ 1000 MeV). At 
these temperatures 6T oc eT'^ldlog^C/dlogT — 1|^^, so 
that the product 7(e, T^cs)ST in Eq. ((37)) is proportional 
to e~^Tj.^J|dlog£/dlogr - 1|^^ Consequently, the con- 
tribution from the low-temperature resonance is dom- 
inant over the contribution from the high-temperature 



resonance. At small values of e (e ^ 0.4), the relic sterile 
neutrino distribution is 



1 



1 — cxp ■ 



dlogC 



dlogT 



(39) 

where k is a constant that depends on sterile neutrino 
parameters and physical constants, and T^es = T^es{^) is 
the low-temperature solution to Eq. ((25)) . This distribu- 
tion can be seen graphically in the lower right panel of 

Fig.m 

Examination of Fig. [1] reveals that Ti.os(e) is a slowly 
decreasing function of e {e.g., a factor of 4 increase in 
e leads to a 15 percent decrease in T^s). As a result, 
€~^T~J remains a decreasing function of e. In addi- 
tion, dlog£/dlogT is maximal immediately preceding 
the resonance cessation event and decreases with increas- 
ing temperature. It follows that |(ilog£/dlogr— 1|~-^ is 
also a decreasing function of e. Thus, the cumulative 

i^e conversion probability, 1 — e"/''"'*^, is a de- 

creasing function of e. The interplay between this and 
the increasing spectrum of Ve produces the peak seen at 
low values of e. 

The high-e peak corresponds to the maximum value of 
e (cmax) for which MSW resonances exist. At this value 
of e, the MSW resonance equation, Eq. ((25)) , has only one 
positive temperature solution (a double root). In the case 
shown in Fig. [1] this occurs at Cmax ~ 1.15. 

For e < eniax, the energy spectrum increases as e in- 
creases, and as it approaches the maximum value, it in- 
creases at a greater rate until it hits the peak. Here, the 
proportionality to e loses the sway it had at lower values 
of e because the fractional changes are an order of magni- 
tude smaller. In this regime, the high-temperature reso- 
nance moves to lower temperatures where it becomes sig- 
nificant in contributing to the Ve Vs conversion proba- 
bility. Additionally, as the two resonances approach each 
other, the peaks tend to blend with each other, enhanc- 
ing the value of Jjdr, culminating with the merger of 
the two peaks. The sharp increase in indicates the 
point at which this phenomenon becomes significant. 

For values of e slightly larger than Cmax, the energy 
spectrum falls off precipitously. This is because neutrinos 
with e > Emax havc no resonances, and thus sterile neu- 
trino production is suppressed compared to the regime 
with resonant production. 

The mid-e peak corresponds to the last resonant value 
of e. At the temperature where MSW resonances cease 
to exist, there is an abrupt transition from resonant to 
nonresonant sterile neutrino production. During the final 
stages of resonant sterile neutrino production, there is a 
large conversion rate of active neutrinos to sterile neutri- 
nos, producing a large lepton number loss rate, leading to 
a smaller resonance width. Immediately after the reso- 
nance cessation event, lepton number is depleted at a far 
lower rate, as only nonresonant conversion remains. This 
leads to the low-temperature side of the resonance being 
much broader than the high-temperature side. As a re- 
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suit, there is an enhancement in the Vs conversion 

probability that is maximized at the last resonant e. This 
effect is significantly smaller than the effects associated 
with the two peaks discussed previously. 

Figures [2] - |4] show the effects of changing each of 
the three free parameters in our calculation. Figure [2] 
presents the effects of changing the initial lepton number 
in each active neutrino species, but leaving the sterile 
neutrino mass and vacuum mixing angle fixed. Figure [3] 
shows the effects of changing the sterile neutrino mass, 
and Figure 0] shows the effects of changing the vacuum 
mixing angle 6. Each figure displays the same four panels 
as discussed for Fig. [1] The results discussed above serve 
as a guide to understanding the trends in these figures. 

In Fig. [21 we have taken nis — 64keV and sin^ 26* = 
10""'^''. We have varied the initial lepton number and 
show three cases, corresponding to L^q — 0.76 x 10^"^, 
1.1 X 10-3, and 1.3 x IQ-^. 

The resonance cessation event occurs at a higher tem- 
perature for lower initial lepton numbers. We can under- 
stand this result using the resonance condition, Eq. (P7)) . 
Since we assume that initially the lepton number is equal 
in each active neutrino species and we only allow Ve-i^s 
neutrino oscillations, the potential lepton number can be 
written as C(T) = 2Leo + 2Le{T). Most of the initial i^e 
lepton number is depleted by v,, — > conversion, so that 
at the cessation event £ « 2ieO- As a result, lower ini- 
tial lepton numbers mean an earlier (higher temperature) 
loss of resonance. 

The maximum resonant e, Emaxi increases with increas- 
ing initial lepton number. As seen in Eq. (pS)) . the upper 
resonance is proportional to £/T^. This implies that 
the resonant scaled momenta can reach higher values for 
higher initial lepton numbers. For e ^ 2.2 there are more 
i^e available to convert into i>s at a given value of e as e 
increases. 

We must perform the same analysis done above in es- 
timating the cumulative Ve — > Vs conversion probabil- 
ity. Preserving the proportionality to the potential lep- 
ton number and sterile neutrino mass and mixing angle, 
we have 



7(e,Ti.cs)(5r cx e ^T^.J 



d\og£ 



dlogT 



1 



£-^m\ sin2 2e'. 



(40) 



where the assumed sterile neutrino mass is m. 



The relevant resonance temperatures are approximately 
equal to the temperature of the resonance cessation 
event, which increases with decreasing initial lepton num- 
ber. In our calculations, T^-J dominates over so 
in concert with the trend for Cmaxi increasing the initial 
lepton number leads to a larger (in magnitude) loss of Vf. 
lepton number, 

^L.^-^^[ U^)d.. (41) 

Consequently, the number density of sterile neutrinos 
produced, ALg, increases with increasing ini- 

tial lepton numbers. Thus, the sterile neutrino closure 



fraction, fJ^ oc rUsn^^ , also increases with increasing ini- 
tial lepton number. 

In Fig. [21 we have taken L^o = 1.1 x 10^"^, 
sin^ 26 = lO""'^" and varied the sterile neutrino rest mass. 
In this figure we show the results of our calculations for 
sterile neutrino masses = 45 keV, 64 keV, and 79 keV. 

It is clear from this figure that the resonance cessation 
event occurs at higher temperatures for higher masses. 
This trend can be understood by noticing that the right 
hand side of Eq. (|?7)) is proportional to nis. Therefore, 
with the same arguments as above, but now with a fixed 
lepton number, equality in Eq. l[27)) (which marks the 
end of resonances) is reached at higher temperatures for 
higher assumed sterile neutrino masses. 

In this calculation, at high temperatures 
(T > 2000 MeV) the resonant values of e do not 
depend on sterile neutrino mass. Using the propor- 
tionality in Eq. (j40p . lepton number depletion, and the 
concomitant sterile neutrino production, will be highest 
for the higher mass sterile neutrinos, at least early on. 
As a result, the highest mass sterile neutrino will have 
the lowest value of Cmax and the highest temperature at 
which the resonance cessation event occurs. 

The two competing factors in Eq. (|^ . T^^J and m^, 
nearly balance each other, but the higher value of Cmax 
for lower sterile neutrino masses results in a larger (in 
magnitude) value of ALg. However, it can be seen in 
the lower left panel of Fig. [3] that the values of ALg for 
the three sterile neutrino masses shown do not differ by 
much. As a result, the sterile neutrino closure fraction, 
the relevant quantity in searching for a dark matter can- 
didate, increases for increasing assumed sterile neutrino 
mass. 

In Fig. [H we have taken nis = 64keV and L^o — 
1.1 X 10"'^ but have employed different values of the vac- 
uum mixing angle 9. The figure shows three cases, for 
sin^ 26 = 0.5 X 10~^", 1.0 x 10^1°, and 2.0 x lO^^". 

For the vacuum mixing angles that we are concerned 
with here, the locations of the MSW resonances are in- 
sensitive to the value of 9 (since cos 26* w 1 — 0.5 sin^ 29). 
As a result, differences in the evolution of e^cs are linked 
to the differences in the lepton number evolution. At 
high temperatures (T ^ 1500 MeV for this calculation), 
the resonances sweep through the Ve distribution at the 
same rate, but the proportionality to sin^ 29 in Eq. (|^D|) 
leads to a larger sterile neutrino production for larger vac- 
uum mixing angles. As a result, the resonance condition, 
Eq. (|27)l . will be violated at the highest temperature for 
the highest mixing angle. Thus, the largest mixing angle 
corresponds to the highest resonance cessation tempera- 
ture and the smallest value of emax- 

The combination of the proportionality of the cumu- 
lative Ve — *■ t^s conversion rate to sin^ 29, along with the 
weak dependence of the resonant temperatures and Cmax 
on the vacuum mixing angle, result in resonant produc- 
tion of sterile neutrinos being highest for higher values of 
vacuum mixing angle. However, this is just the tip of the 
iceberg. Nonresonant scattering-induced decoherent pro- 
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FIG. 2: The four panels present the same quantities as in Fig. [T] In this figure, we have taken nis = 64keV, sin^ 26 = 10^^", 
and three different values for the initial lepton number: L^o — 0.76 x 10"'' (short dashed curves), 1.1 x 10"'^ (long dashed 
curves), and 1.3 x 10~^ (solid curves). 



duction of sterile neutrinos is also proportional to sin^ 29. 
As a result, higher vacuum mixing angles generally result 
in a higher sterile neutrino closure fraction. 

In Fig. [51 we adopt the same parameters as in Fig. [1] 
{ms = 64keV, sin^ 26* 10"^°, and Leo = 1.1 x 10"^), 
but we modify the calculation by increasing the total 
scattering rate, Fg. In the calculations for Figs. [1] - 21 
we have used a value of j/e = 1.27 in the total scattering 
rate, Eq. (fTB)) . However, as discussed in Sec. Ill Al this 
is a value that is relevant for temperatures much lower 
than those in our calculations. 

We modify our calculations by assuming that the to- 
tal scattering rate is proportional to the number of weak 
degrees of freedom that are in thermal equilibrium in the 
early universe. The value y,, « 1.27 stems from Ve scat- 
tering on pairs and neutrinos and antineutrinos of all 
three active flavors. However, for the temperatures that 
we are concerned with here, there are also fi^ pairs and 
up, down, and strange quarks and antiquarks thermally 
populated in the early universe. By including all three 
quark colors and the fact that the weak current is left- 
handed, we modify the total scattering rate by a factor 
of 3.5. In other words, in testing the sensitivity to the 
scattering rate, we substitute Te{e,T) —^ 3.5re(e,r). 

This calculation shows that there is little effect in the 
resonant conversion regime. This result can be under- 
stood by noticing that 7(e,T',„ax) oa and ST oc Fg. 



Therefore, the resonant production of sterile neutrinos is 
only weakly dependent on the total scattering rate. In 
Fig. [5] this is readily apparent, because the evolution of 
the resonant scaled momentum and the lepton number 
(and consequently the sterile neutrino production his- 
tory) are nearly identical for T ^ 950 MeV. 

The most dramatic effect of an increased scattering 
rate can be seen in the sterile neutrino production his- 
tory after the resonance cessation event. Roughly 55per- 
cent of the final sterile neutrino relic density is produced 
nonresonantly (compared with 30 percent in the lower Fe 
case) . This is because away from resonance the quantum 
Zeno effect is less efficient at suppressing sterile neutrino 
production, allowing the Ve ^ rate to increase with 

Another effect of the increased scattering rate is seen 
in the sterile neutrino energy distribution function. The 
most significant difference is that the high-e peak is not 
as sharply peaked in this case. This is due to the com- 
bination of two processes. Nonresonant sterile neutrino 
production for e > Cmax is enhanced with a larger scat- 
tering rate, producing a less steep drop off. The merging 
of the two resonances that created the sharp peak in the 
low Fe case is less effective for a larger Fg. As discussed 
above, the peaks are broader with a smaller amplitude. 
Consequently, when these peaks come together, the re- 
sult is a smaller cumulative sterile neutrino production. 
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FIG. 3: The four panels present the same quantities as in Fig.[T] In this figure, we have taken Leo ~ 1-1 x 10"^, sin^ 29 — 10^^", 
and three different values for the sterile neutrino rest mass: = 45keV (solid curves), 64keV (long dashed curves), and 79keV 
(short dashed curves). 



To this point we have only discussed sterile neutrino 
dark matter candidates that are ruled out by the cur- 
rent X-ray observations. In Fig. [6] we present the energy 
spectrum for a sterile neutrino dark matter candidate 
that is still "in play" as far as the x-ray observations are 
concerned. To evade the x-ray limits, we must reduce the 
sterile neutrino mass and/or vacuum mixing angle. From 
the studies we conducted above, we know that the relic 
sterile neutrino closure fraction is more strongly depen- 
dent on mass than on mixing angle. Therefore, we must 
reduce both mass and mixing angle to get under the x- 
ray limits. However, to get the correct relic density we 
also increase the initial lepton number and the scattering 
rate as discussed above for Fig. [5] 

For this calculation we use a sterile neutrino rest 
mass TOs — 12keV with vacuum mixing angle sin^ 26 = 
3.2 X 10~^^ and initial lepton number LeO = 4.8 x 10~^. 
The sterile neutrino properties are just below the cur- 
rent x-ray constraints, while the initial lepton number is 
much less than the best constraints on the lepton num- 
bers in the early universe. We display only the energy 
spectrum because the other three panels do not provide 
much physical insight in this case. 

For this case our calculations give e^ax ~ 200, which 
means that the resonance sweeps through the entire 
distribution. As a result, the high-e peak cannot be seen 
in Fig. [HI The mid-e peak is at e « 60, so that is also 



missing. We are left with the low-e peak and a long 
tail at higher e that eventually disappears because of the 
scarcity of neutrinos in the v,, distribution at these scaled 
momenta. 

For this particular case the resonance cessation event 
occurs at T « 100 MeV, but as discussed above, this cor- 
responds to e ~ 60. The low-temperature resonance for 
more relevant values of e (i.e., e ~ 3) occurs at tem- 
peratures 170 MeV. Since the resonances persist to 
such low temperatures, there is very little nonresonant 
production of sterile neutrinos. As a result, the relic 
sterile neutrino density produced is relatively insensitive 
to the total scattering rate. However, this temperature 
is of interest in the early universe as it is the epoch of 
the QCD transition. The nature of this transition affects 
the specifics of the sterile neutrino production [i^ but 
in broad brush should serve to enhance the final sterile 
neutrino closure fraction. 

Also plotted on Fig.[n]is a Fermi-Dirac spectrum whose 
normalization is chosen to have an equivalent sterile neu- 
trino density. The peak in the sterile neutrino spectrum 
is at e « 1, while the Fermi-Dirac spectrum peaks at 
e « 2.2. We see that a significant fraction of the total 
sterile neutrino density has smaller values of e than in its 
Fermi-Dirac counterpart. Nonetheless, the high-e tail of 
the distribution yields an average value of e, (e) « 2, com- 
pared to an average e of 3.15 for a Fermi-Dirac spectrum. 
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FIG. 4: The four panels present the same quantities as in Fig. [T] In this figure, we have taken rris = 45keV, Leo = 1.1 x 10 , 
and three different values for the vacuum mixing angle; sin 261 = 0.5 X 10""' (solid curves ), 1.0 X 10"^° (long dashed curves 1 , 
and 2.0 x 10'^° (short dashed curves). 



As a result, the sterile neutrino population produced in 
our calculations is "colder" by a factor of about 1.6 than 
an equivalent population with a Fermi-Dirac spectrum. 
This is qualitatively consistent with the findings in Ref. 
[l6| . but the sterile neutrino population here is not as 
cold. 

Our numerical calculations show that the sterile neu- 
trino rest mass, vacuum mixing angle, and initial Icp- 
ton number parameters which produce the correct dark 
matter closure fraction, but which evade X-ray bounds, 
produce relic Vs energy spectra which are qualitatively 
similar to that in Fig. [S] Since this generic spectrum 
is skewed relative to a Fermi-Dirac form toward lower 
values of e, viable dark matter sterile neutrinos which 
are produced via lepton number enhancement should be 
"colder" than their rest masses would suggest. 

Finally, in Fig. [7] we perform a consistency check to 
determine the validity of the quantum Zeno ansatz that 
we have used throughout our calculations. To do this, 
we use the sterile neutrino parameters and initial lepton 
numbers employed in the calculations for Fig. [1] (m^ = 
64keV, sin^ 261 = 10"^", and L^o = 1.1 x 10"^), and also 
use the smaller total scattering rate discussed above. We 
solve the quantum kinetic equations presented in Sec. 
Ill Al for e = 1, using the lepton number evolution from 
our calculations. 

Plotted in Fig. [7] are the instantaneous — i^s con- 



version rate per unit temperature interval, \dfs/dT\ (solid 
curve), and the fractional difference between the rate in- 
ferred from our solution to the quantum kinetic equations 
and the rate given by the quantum Zeno ansatz (dashed 
curve). Note that the actual conversion rate is always 
slightly larger than the rate from the quantum Zeno ef- 
fect, but in general the two agree within a few percent. 
Note also that the differences change by an order of mag- 
nitude at the resonances. 



Since the actual rate is possibly slightly higher than the 
rate we used with the quantum Zeno ansatz, we probably 
under-produce sterile neutrinos in our calculations. The 
interval between the two resonances is an important one 
in determining the final sterile neutrino production. If 
corrections on the order of a percent occurred for all our 
values of e, this could increase the temperature of the res- 
onance cessation event and, in turn, increase the ultimate 
relic sterile neutrino density, possibly in a nonlinear fash- 
ion. Our numerical calculations show that a 1 % increase 
in the conversion rate gives a 0.3 % increase in relic ster- 
ile neutrino density. We conclude that a small change in 
the conversion rate does not lead to a significant change 
in the relic sterile neutrino density. 
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FIG. 5: The four panels present the same quantities as in Fig. [T] In this figure, we have taken rris = 64keV, Leo ~ 1-1 x 10"'^, 
and sin^ 20 = 10"^*^. We present two calculations, one with a low total scattering rate Fe as discussed in section fll Al (solid 
curves), and one with a higher total scattering rate Fe — > 3.5Fe (dashed curves). For T > 950 MeV, the cumulative sterile 
neutrino production history, lepton number history, and resonant values of e are nearly identical in the two calculations. 
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FIG. 6: The sterile neutrino energy spectrum as a function of 
scaled momentum e for an "in play" sterile neutrino dark mat- 
ter candidate with rus = 12keV, sin 29 = 3.2 X 10"", and 
Leo — 4.8 X 10^'^ (solid curve). Also shown is a Fermi-Dirac 
spectrum normalized so that both spectra would produce the 
same number density of sterile neutrinos (dashed curve). 
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FIG. 7: The solid curve is the sterile neutrino conversion 
rate per unit temperature interval, \dfs/dT\ for e = 1 in 
our calculations with nis = 64keV, sin 261 = 10"^^ and 
Leo = 1.1 X 10"^. The dashed curve shows the fractional 
difference between the rate inferred from a solution of the 
quantum kinetic equations to the rate implied by the quan- 
tum Zeno ansatz. 



V. CONCLUSIONS 

One option for producing viable sterile neutrino dark 
matter is lepton number enhancement of scattering- 
induced decoherence. The calculations presented here 
address the basic physics of this production process. A 



key conclusion of our work is that the interplay of lep- 
ton number depletion and resonant sterile neutrino pro- 
duction leads to a set of generic peak features in the 
relic sterile neutrino energy spectrum. We have stud- 
ied how these peaks depend on the assumed initial lep- 
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ton number, the magnitude of the assumed active-sterile 
couphng (i.e., the vacuum mixing angle), and the sterile 
neutrino rest mass. We conclude that relic sterile neu- 
trinos with parameters that evade x-ray constraints will 
have a characteristic single-peak energy spectrum. This 
spectrum will be "cold" compared to a Fermi-Dirac form 
spectrum with the same integrated relic number density. 
This could relax to some extent Lyman-alpha forest con- 
straints on sterile neutrino rest mass [i^l . 

We have compared the widely-used quantum Zeno 
ansatz with a full quantum kinetic equation treatment 
of scattering-induced decoherence production of sterile 
neutrinos in the early universe. Our conclusion is that 
the quantum Zeno approximation is adequate, at least 
during the epoch in the early universe when ~ keV mass 
sterile neutrinos are produced. 

Finally, we should note the limitations of this study 
and point out the topics that need further detailed ex- 
amination. The calculations presented here are crafted 
to study the physical effects which determine relic sterile 
neutrino energy spectra. They only treat schematically 
the microphysics of scattering degrees of freedom and fi- 
nite temperature in-medium effects. We have attempted 
to gauge the sensitivity of final relic Vg densities and en- 
ergy spectra to these issues by simply varying the pre- 
scription for active neutrino scattering rates. Though we 
find that our generic spectral features are relatively in- 
sensitive to changes in these rates, our calculations never- 
theless cannot serve to predict relic densities accurately. 

At least three aspects of our calculations for sterile 
neutrino production and lepton number depletion are 
simplistic: (1) we have allowed only a single channel, 
Ue i^s, for sterile neutrino production; (2) we have 
taken initial, preexisting lepton numbers in all active neu- 
trino species to be the same; and (3) we have not allowed 
mixing between active neutrino flavors. As for the last 
point, the study in Ref. [l^ suggests that dynamically 



including active-active mixing should be roughly similar 
to having a larger initial value of LeO in our model. The 
idea being that as are converted to Ug, some j/^ and 
Vt neutrinos convert into v,,- In essence then, the net 
lepton number in the mu and tau neutrino seas serves 
as a reservoir that "feeds" the sea as sterile neutrinos 
are produced. In fact, the measured active-active mixing 
parameters may suggest mixing between the flavors is ef- 
ficient and that equality of lepton numbers is a plausible 
assumption [2^ . A larger reservoir of active neutrinos 
feeding into Uf. Vg would not change the qualitative as- 
pects of our conclusions on the interplay of lepton num- 
ber depletion and neutrino asymmetry-enhanced sterile 
neutrino production in the early universe. 

In contrast, if there were more than one kind of ster- 
ile neutrino, or if there were multiple channels of active- 
sterile conversion {e.g., v^j^^r ^ i^s m addition to i^e i^s), 
then both lepton number depletion and sterile neutrino 
production histories might be different than those cal- 
culated here. In fact, a key conclusion of our work is 
that these histories, as well as the relic sterile neutrino 
energy spectrum that goes with them, are a product of 
the potentially complicated coupling between the sterile 
neutrino production channels and the way lepton num- 
bers evolve. This is a rich problem which deserves further 
consideration. 
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